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Abstract 

The quantum electrodynamical (QED) short wavelength correction on plasma wave 
propagation for a non-relativistic quantum plasma is investigated. A general dispersion 
relation for a thermal multi-component quantum plasma is derived. It is found that the 
classical dispersion relation for any wave mode can be modified to include quantum 
and short wavelength QED effects by simple substitutions of the thermal velocity 
and the plasma frequency. Furthermore, the dispersion relation has been modified 
to include QED effects of strong magnetic fields. It is found that strong magnetic 
fields together with the short wavelength QED correction will induce dispersion both 
in vacuum and in otherwise non-dispersive plasma modes. Applications to laboratory 
and astrophysical systems are discussed. 

1 Introduction 

There is currently a great deal of interest in quantum vacuum effects (see e.g. Refs. [1-10] 
and references therein for an up-to-date selection) as well as in non-relativistic quantum 
plasma systems [11-23], starting from e.g. the Schrodinger model of the electron [24,25] or 
the Pauli equation [26,27]. Using a theory based on non-relativistic quantum mechanics, 
there have been studies of a wide range of plasma wave phenomena [13, 15, 16, 18] as well 
as astrophysical applications (e.g. Ref. [17]). In relativistic quantum theory, some of 
the literature has been motivated by the rapid progress in experimental techniques [28-30] 
suggesting experiments for different vacuum polarization effects, [2-4,7,31-38] (see also Ref. 
[9] for a discussion). For example, in Refs. [2,7,31,32,34-38] suggestions for experimental 
setups to verify photon-photon interaction are given, including the use of strong lasers (e.g. 
Ref. [2]) and dense plasma columns [36]. New experimental results have been presented 
in for example Ref. [3], which reports a change in polarization of light traveling through 
a strong magnetic field (see also Refs. [4,10] for further discussions). Other examples are 
Ref. [39] which reports electron-positron pair production due to a intense electromagnetic 
field interacting with high frequency photons, and Ref. [40] which investigates non-linear 
Compton scattering. There are also studies of theoretical and/or observational nature 
[5,6,41-54], which find applications in high intensity laser systems [5,6,51,52] and in 
astrophysics [41-50,53,54]. For quantum plasma dynamics, a set of effective macroscopic 
field equations may be derived from the microscopic equations of motion for the relevant 
particles using a range of techniques [24,25,28]. Such effective models can be very useful 
when trying to understand some of the many collective effects that may occur in such 
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systems, and can also give an indication on when one can expect macroscopic phenomena to 
arise from the quantum realm. Astrophysical environments exhibit many extremes [55,56], 
e.g. magnetars [57], and as such constitute interesting areas of applications for quantum 
vacuum physics as well as relativistic quantum plasmas. For such applications there have 
been a number of studies of relativistic quantum plasmas, e.g. Refs. [58-61]. Indeed, in 
order to properly investigate such extreme systems, a quantum field theoretical basis that 
incorporates certain quantum vacuum effects is desirable. 

Here, we will investigate the effects on non-relativistic quantum plasmas due to short 
wavelength corrections from quantum electrodynamics (QED). Such combined effects are 
indeed compatible, as the short wavelength corrections do not necessarily imply relativis- 
tic particle motions. In Section 2 we present the weak field short wavelength corrected 
Maxwell's equations and derive a quantum corrected fluid model of the plasma. These 
modified equations are used in Section 3 to derive a general dispersion relation for any 
plasma mode in a thermal multi-component plasma. The possibility of designing experi- 
ments where some of the above effects could be detected is of great fundamental interest. 
Such scenarios are discussed in Section 4. Furthermore, the quantum field theoretical cor- 
rection due to strong magnetic fields, relevant to astrophysical environments, is discussed 
in Section 5 where the general dispersion relation is modified to also include such effects. 
In particular, we derive the dispersion relations for the possible plasma modes in environ- 
ments similar to those in the vicinity of magnetars and pulsars. Finally, in Section 6 we 
summarize our results. 



2 Governing equations 

2.1 Electromagnetic field equations 

The first order QED effects can effectively be modeled through the Heisenberg-Euler La- 
grangian density [62,63]. This Lagrangian describes a vacuum perturbed by a slowly 
varying electromagnetic field. The effect of rapidly varying fields can be accounted for by 
adding a derivative correction to the Lagrangian [64]. This correction is referred to as the 
derivative QED correction or the short wavelength QED correction. Since the QED effect 
of electron-positron pair creation is not included in this model, the field strength must be 
lower than the critical field strength, E , cr ; t = m^c 3 /he « 10 16 V/cm -1 , and the frequency 
must remain lower than the Compton frequency to e = m e c 2 /H [28,65]. Here m e is the 
electron mass, c is the speed of light in vacuum, e is the elementary charge and H is Plank's 
constant. The Heisenberg-Euler Lagrangian density with the derivative correction reads 

C = Co + Che + Co 

= £ j-F ab F ab + £ ^r 4(F ab F ab Y + 7[F ab F ab y + ae 
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where Co is the classical Lagrangian density, while Che represents the Heisenberg-Euler 
correction due to first order strong field QED effects, Co is the derivative correction, 
□ = d a d a is the d'Alembertian, F ab is the electromagnetic field tensor and F ab = e abcd F ca -/2 
where e abcd is the totally antisymmetric tensor. The parameter k = 2a 2 ^, 3 /45m 4 c 5 gives 
the nonlinear coupling, a = (2/15)ac 2 ju\ is the coefficient of the derivative correction and 
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a = e 2 /Aithceo is the fine structure constant, where £o is the free space permittivity. We 
obtain the field equations from the Euler-Lagrange equations df, [d£/dF a t,] = poj a [66,67], 



{\ + 2aU)d a F ab = 2e Kd a 



F pod ) F ™ + ±i Fcd F ca ) F 



rid,/) 



+ W)3 



(2) 



where j a is the four-current and /j,q is the free space permeability. 

For the remainder of this manuscript we will investigate low-amplitude linear wave 
perturbations in a magnetized quantum plasma. The strong field QED term proportional to 
k (F c dF cd ^ F ab then contributes with a term that is quadratic in the unperturbed magnetic 
field .Bo, and linear in the wave perturbation. However, due to the smallness of k such 
a term does not become dynamically important unless the unperturbed field strengths 
approaches Bq ~ 10 8 — 10 10 T, i.e. in pulsar or magnetar environments. An example of 
such effects will be considered in the end of section 5, but until then the strong field QED 
effects will be omitted. Following Ref. [68] we will keep the derivative QED corrections 
proportional to a. This is a natural approach provided cBQ/E crit -C u/uj e , where w is a 
typical frequency of the field. The corresponding sourced Maxwell equations resulting from 
the derivative corrected field equation then become 



i d 2 



V-E = ^, 



1 + 2a 



c 2 dt 2 
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and the source free Maxwell equations are V • B = and 

V x E 



dB 



(3) 



(4) 



(5) 



Here j = q s n s u s is the current density, p = q s n s is the charge density, the index s 
denotes the species of particles considered, q s is their charge, n s is the particle density and 
u s is the particle velocity. 



2.2 Particle equations of motion 

The plasma can be described by a model which captures some of the quantum properties 
of the plasma particles. We here follow Ref. [24], where the particles are described by the 
statistical mixture of N states V^j i = 1,2, ...,iV. The index i sums over all particles 
independent of species. The one-particle states satisfy the Schrodinger equation, which for 
each state reads 

= ~2~m^ 2 ^ i + qi ^ u ( 6 ) 

where rrii and qi are the mass and the charge of the particle respectively. The equations 
are coupled by Poisson's equation 

1 N 

vv= — y>^i 2 (7) 

where pi is the occupation probability of state ipi. This model amounts to assume that all 
entanglement between particles are neglected. To derive a fluid description we make the 
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ansatz ipi = ^/riiexp(iSi/h) where re, is the particle density, Si is real, and the velocity 
of the i'th particle is Uj = VSi/nii. Next we define the global density and velocity for 
particles of species s as n s = J2jPj n j an d u s = Ylj Pj n j u j/ n s, where j runs over all 
particles of species s. Inserting this into equations J6]) and and separating the real and 
the imaginary part, we obtain the continuity equation 

+ V • (n s u s ) = 0, (8) 

and the momentum equation 

^ + K • v) u s = — e — -J-vp + * v (J-v*^r\ . (9 ) 

at m s m s n s lm A s \\Jn s J 

Here we have made the substitution J2jPjC^ 2 \/^j)/\/^j ~^ ( v v^sVv^s) which is valid 

1 /3 

for length scales larger than the Fermi length Xf s = vf s /uj ps , where vf s = fi {3Tr 2 £ou 2 s /m 2 q 2 ^ 
and ujp S = \J q^nQ S /eom s is the Fermi velocity and the plasma frequency of species s. We 
stress that the pressure term contains both the fermion pressure, Pp S; and the thermal 
pressure, Pf. 

This model can be extended to include non-zero magnetic fields, by starting with a 
Hamiltonian which includes the vector potential. In the resulting model [14], the continuity 
equation ([8]) is unchanged, and the momentum equation (JU) becomes 

' + (u s • V) u s = — (E + u,xB) - —VP + -^V ( -Lv 2 ^;) . (10) 



dt m s m s n s 2m 2 . \^/n, 

3 The general dispersion relation 

Our next aim is to study the linear modes of the system described by eqs. (SJ) , (EJ) , (ED and 
(fTUI) . We linearize and Fourier decompose these equations and let tiq s represent a constant 
unperturbed particle density while n\ s denotes a density perturbation. Furthermore, the 
unperturbed velocity is put to zero and from now on u s denotes the wave perturbed 
velocity. For long wavelengths, such that hk/m <C vf *C u/k, the pressure term in eq. 
((EJ) can be written as (l/m s n s )V (P ts + P Fs ) = (l/no«) {v 2 s + | v Fs ) Vn is [24], where 
we have neglected the thermal pressure that particles of different species exert on each 
other. Here, k is the perturbation wave number, u is the perturbation frequency and vt is 
the thermal velocity. The pressure term and quantum correction in eq. (fTOl) can now be 
collected into an effective term with an effective wavenumber dependent velocity V^, 

Below, the collective thermal and quantum effects described by V will sometimes be re- 
ferred to as effective thermal effects. We define the z-axis to be along Bo and the y-axis to 
be orthogonal to k, such that Bo = Bqz and k = k x x + k z z = k±x + k\\z. Now use eq. (JHJ) 
to express n± s in terms of u s , and insert this result into eq. (fTUI) . Next we rewrite eq. (fTOl) 
in the form cr" 1 • u s = E, and find <j s such that u s = a s • E. Using this result together 
with eq. j5|), we can then rewrite eq. (Jl| in the form D • E = where 
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(12b) 



Here I is the identity matrix, ^ s % s is the short wavelength QED and quantum corrected 
plasma susceptiblitiy, uj cs = q s Bo/m s is the gyrofrequency, 



C = 2a [-=■ - k 



2au) 2 , 2 x 



(13) 



gives the derivative QED correction due to rapidly varying fields, n± = k±c/ui and nu = 
k\\c/uj is the index of refraction orthogonal respectively parallel to the magnetic field such 

that n = ^ In 2 , + n?, is the total index of refraction, and we define 
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k 2 V 2 ] 



k 2 V 2 



k 2 , k 2 Vf. 



(14) 



Setting the determinant of D equal to zero gives the dispersion relation for all plasma 
modes. From expression (|12bp it becomes transparent that the dispersion relation for any 
classical plasma mode can be modified to include the quantum correction by making the 
substitution 



Vts 



fl5a) 



and the short wavelength QED contribution can be included by making the substitution 



CO 



ps 



i-C 



(15b) 



Note that the substitution of u; ps should not be made within the expression for vp s . 

In order to demonstrate the usefulness of our results (|15al) and (I15b|) and verify their ap- 
plicability, we compare with previous results. As is well known, the classical dispersion rela- 
tion for ion-acoustic waves with zero ion-temperature reads u/k = vt e /\/l + k 2 v 2 e m e /m,i, 
where indices e and i denotes electron and ion, respectively. Making the substitution 
(|15ap in the classical expression immediately give the dispersion relation (23) of Ref. [13]. 
Furthermore, the classical dispersion relation for almost perpendicular propagating electro- 
static ion-cyclotron waves reads J 2 = uj 2 ci + k\v 2 e m e /mi. Again, making the substitution 
Qloap in the classical formula we obtain agreement with the dispersion relation (13) in 
Ref. [22] in the appropriate limit. 

Moreover, the classical zero temperature dispersion relations for plasma oscillations 
a; 2 , electromagnetic waves in an unmagnetized plasma uJ 2 = k 2 c 2 + lj 2 electromag- 



netic waves propagating parallel to an external magnetic field (u±u) c )(oj—k c /u 



and extraordinary waves (u 2 — uj 2 ){uj 2 + k 2 c 2 — uj 2 ) + uj 2 {k 2 c 2 — uj 2 ) = 0, all give agree- 
ment with corresponding relations in Ref. [68], when substitution (|15bj) in the classical 
expressions. 
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4 Laboratory application 



It would be of great fundamental interest if the quantum or derivative QED effects could 
be detected in experiments. Firstly, we will be concerned with the possibility of detecting 
the derivative QED contribution in the laboratory. To single out the small dispersive QED 
effects from the general dispersive effects of a plasma requires detailed knowledge of the 
plasma parameters for an experiment to be conclusive. We now investigate whether it 
may be possible to detect the short wavelength QED effects by studying how the index 
of refraction scales with frequency close to a plasma resonance, rather than trying to 
determine the actual phase shift caused by the derivative QED correction. For this purpose, 
assume that we let several laser beams with slightly different frequencies pass close to each 
other through the same plasma column in such a way that all beams experience the same 
plasma density perturbations. In this way, precise knowledge of the plasma density profile 
is not necessary for interpreting the results. The phase shift due to the difference in optical 
path length between a beam passing through the plasma column and one that does not 
pass through the plasma column can be measured with high precision [69]. How the index 
of refraction scales with the frequency can be estimated by comparing the phase shifts of 
the different laser beams. 

Since the short wavelength QED effect is enhanced for large wavenumbers, see Ref. [68], 
one of the most interesting modes when considering the derivative QED correction is the 
extraordinary mode (X-wave). The classical dispersion relation for an X-wave in a thermal 
electron plasma reads 



2 1 _ fp_ v vi 

UJ 2 (uj 2 - k 2 v 2 -uj 2 - w p) ' 

We include the quantum and the derivative QED contributions, by making the substitu- 
tions given in (|15ap and (|15b|) . The index of refraction to the lowest non vanishing order 
in a and V then becomes 

^ 1 t ^ y , u ( (uj 2 -1)(uj 2 -2) (uj 2 -lf \ 



2 (uj 2 - k 2 v 2 - ul) 




where the parameters have been normalized according to 

- - 2uJ p t7 - V - - u a - - UJc ns^ 

a = — ^-a, V — — , u = — and uj c = — . (18) 

C Z C LOp LOp 

Here we have disregarded the /c-dependence of the effective thermal velocity and treated 
V as a constant. The classical dispersion diagram for an X-wave has been illustrated in 
Fig. JT]). The derivative QED contribution and the effective thermal contribution given in 
eq. (fTT|) has been illustrated separately with a and V normalized to one. 

For an X-wave close to plasma resonance, eq. (fT7l) shows that the leading term in 
the derivative QED contribution to the index of refraction has a cubic dependence on the 
resonance term, (uj 2 — Oj 2 — l) , in contrast to the classical contribution which is linear 
in this term. For detection to be conclusive, the thermal effects must be much smaller than 
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(c) 

Figure 1: The derivative QED correction (dashed line) and the thermal and quantum 
correction (solid line) to the index of refraction, n, for an X-wave (dash-dot line). Here, 
a and V are both normalized to one and the magnetic field strength is chosen so that 
u) c = 0.5 in Fig. la), uj c = 0.2 in Fig. 16) and u c = 0.05 in Fig. lc) 
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the derivative QED effects in the regime of detection. Unfortunately, the leading term for 
the thermal contribution has the same form as that of the derivative QED contribution. 
This gives us the condition V < uj p \f2a, which must be fulfilled for the thermal effects to 
be smaller than short wavelength QED effects close to a plasma resonance. Even if the 
plasma temperature is low, the Fermi velocity (oc u> p ) will always remain greater than 
u> p V2a for laboratory plasma densities. Thus, derivative QED effects will be negligible 
compared to the effective thermal effects in all laboratory plasmas. 

Even if we allow for higher densities, we can still expect the short wavelength QED con- 
tribution to be negligible. Moreover, care has to be taken since high plasma densities may 
imply a relativistic Fermi velocity, while we in our model have assumed a non-relativistic 
plasma. In contrast to being close to resonance, we next consider a regime far from res- 
onance. In particular, we require Q <C 1. For to to be greater than the cutoff frequency 
while the condition uj -C 1 is satisfied, we must require u) c 3> 1. The leading term of 
the thermal contribution and the derivative QED correction in (fT7l) approximately be- 
comes V 2 /id^id 2 and 4ct/u) 4 u; 2 , respectively. Hence, the condition under which the short 
wavelength QED contribution is greater than the effective thermal contribution becomes 
approximately V < uJ p V^a. Thus, derivative QED effects will be negligible compared to 
effective thermal effects also in this regime. 

One could in principle try to use the same technique outlined above to detect the 
dispersive effects due to the quantum correction. With the quantum correction included, 
the effective velocity V 2 , should be written as 

y2 = v 2 +n 2 (19) 

where v 2 = ^v 2 F + v 2 . In order to get a phase shift due to the quantum contribution of 
approximately it in a 10 m long plasma column with a plasma frequency in the infrared 
regime (u p as 10 14 Hz), v = v/c must be smaller than 10 -7 for the quantum correction 
to be the leading effect. The Fermi velocity alone gives a much greater contribution at 
these densities, vf/c ps 10~ 4 . Furthermore, the thermal effect removes the resonance 
in the dispersion diagram at high temperatures. In this regime, the quantum correction 
consequently stretches across the resonance as can be seen in Fig. . 



5 Strong field QED effects 

Another interesting regime is that of extreme magnetic fields. In the vicinity of magnetars 
and pulsars, the magnetic field Bq can reach 10 10 — 10 11 T close to the surface. Returning 
to the full wave equation (j2J) accounting for strong field effects, we note that the term 
proportional to k contributes with terms that are linear in the wave field and quadratic 
in Bq. Following Ref. [70] we find that the tensor in expression (I12a|) can be modified to 
include the QED effect of strong magnetic fields by adding the correction 

c / 1 — ra| n±n\\ \ 

Xqed = - 4 T37 l~n 2 -2n L , (20) 

V "±«|| ~\-n\ ) 

where £ = k,£qc 2 B 2 = (a/907r)(c.Bo/i?crit) 2 - In magnetar environments with extreme 
magnetic fields we expect the presence of a electron-positron plasma [55]. Assuming n e o ~ 
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Figure 2: The quantum correction to the index of refraction for an X-wave at different 
values of v. The plasma frequency is in the infrared regime (u p = 10 14 Hz) and the 
magnetic field strength is 1 Tesla (p c = 0.001). 



Upo as well as u p , u <C u c , the (1 : 2), (2 : 1),(2 : 3) and (3 : 2) elements of expression (|12bj) 
are small and can be neglected. The dispersion relation for a plasma in a strong magnetic 
field where strong field QED and short wavelength QED corrections are included can now 
be derived from expression (|12a|) with expression (|20l) added. We find that the dispersion 
relation separates into two modes, 



u>l (lo 2 - k 2 V 2 ) 

1 - n 2 (1 - C - 4£) - 8£ni + L „ (21a) 

uj 2 L 2 -k 2 v 2 ^ 



1 I (!-"«) 

u 2 -k 2 V 2 




and 



where to 2 = oj 2 e +u> 2 p . Terms proportional to ui 2 /ui 2 in the second mode has been dropped, 
and £ >C 1 is assumed. Eqs. (I21al) and (|21bj) generalized the dispersion relations derived in 
Ref. [70] to account for quantum and thermal effects, as well as QED derivative corrections. 
We note that if we let ( — ► 0, the group velocity of the mode (|21a|) becomes a constant 
in the zero temperature limit, whereas a nonzero C implies dispersive propagation. While 
the dispersive effects are small for most part of the spectrum, we note that for frequen- 
cies approaching the Compton frequency uj e together with field strengths approaching the 
critical field £7 c rit> dispersion may be appreciable. Strictly speaking, radiation with u > uj e 
propagating in external fields with cBq > E cr n lies outside the validity of our model, and 
the dispersion relations (|21aj) and (I21b|) may thus not be entirely accurate in this regime. 



Nevertheless, our results here show convincingly that such electromagnetic wave propaga- 
tion must be significantly dispersive. For high intensities, radiation in pulsar and magnetar 
environments with the above properties give raise to effects such as pair production and 
photon splitting (see e.g. Ref. [61] and references therein). The competition between these 
two types of processes may determine the pair-plasma density [71] in the vicinity of the pul- 
sars and magnetars, and thereby have observable consequences. Normally, such effects are 
described without accounting for the dispersive wave properties [41,54], but we conclude 
here that our above results indicate that wave dispersion can play a role in this context. 

6 Summary 

In this paper we have investigated how the derivative QED correction and the quantum 
correction affects plasma wave propagation in a multi-component thermal plasma. We have 
derived the generalized dispersion relation for all plasma modes. Specifically, it is found 
that the dispersion relation for any classical plasma mode can be modified to include the 
quantum contribution and the derivative QED correction by simple substitutions of the 
thermal velocity and the plasma frequency according to (|15a[) and (I15b|) . Furthermore, the 
QED effect of a strong magnetic field has also been included in our plasma description. 

We have also investigated the difficulties with detection of the short wavelength QED 
correction and the quantum correction. It is found that dispersive effects due to the Fermi 
velocity in plasma wave propagation will dominate in laboratory regimes. We conclude that 
the detection of the short wavelength QED correction or the quantum correction to the 
dispersion relation in laboratory plasma regimes is difficult. However, there may perhaps 
be ways around these particular problems. 

Finally, we point out that there may exist astrophysical environments where strong field 
QED effects can be of importance. In particular, magnetars and pulsars offer environments 
with extreme magnetic fields. We have derived the dispersion relations for the two plasma 
modes likely to dominate in these environments. The short wavelength QED effect will 
give rise to dispersion in wavemode (|21a)) which otherwise would be dispersionless in the 
zero temperature limit. Furthermore, strong magnetic fields will induce dispersive effects in 
vacuum. This is a unique property arising from the short wavelength QED correction. The 
dispersion can be significant for field strengths close to or above the critical field strength. 
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